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ABSTRACT 



Studies of jets from young stellar objects (YSO's) suggest that material is launched from 
a small central region at wide opening angles and collimated by an interaction with the sur- 
rounding environment. Using time-dependent, numerical magnetohydrodynamic simulations, 
we follow the detailed launching of a central wind via the coupling of a stellar dipole field 
to the inner edge of an accretion disc. Our method employs a series of nested computational 
grids, which allows the simulations to follow the central wind out to scales of tens of AU, 
where it may interact with its surroundings. The coupling between the stellar magnetosphere 
and disc inner edge has been known to produce an outflow containing both a highly colli- 
mated jet plus a wide-angle flow. The jet and wide-angle wind flow at roughly the same speed 
(100-200 km s _1 ), and most of the energy and mass is carried off at relatively wide angles. 
We show that the addition of a weak disc-associated field (<C 0.1 Gauss) has little effect on 
the wind launching, but it collimates the entire flow (jet + wide wind) at a distance of several 
AU. The collimation is inevitable, regardless of the relative polarity of the disc field and stel- 
lar dipole, and the result is a more powerful and physically broader collimated flow than from 
the star-disc interaction alone. Within the collimation region, the morphology of the large- 
scale flow resembles a pitchfork, in projection. We compare these results with observations of 
outflows from YSO's and discuss the possible origin of the disc-associated field. 
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1 INTRODUCTION 

Several decades of observing young stellar objects (YSO's) 
have revealed a ubiquitous coincidence of structured outflows 
and accre tion discs, accompanyin g the star for mation pro- 
cess. See | Reipurth & Baihl i200ll) . lEisloffel et all <200Ch . and 
iKQnigl & PudritzH200d) for detailed reviews of Herbig-Haro (HH) 
objects and YSO outflows in general. The outflows are invari- 
ably collimated to some degree and always travel nearly perpen- 
dicular to accretion discs. Often, relatively wide-angle molecular 
(e.g., CO) outflows escort narrower optical outflows from the same 
source. The optical outflows typically consist of jets, which are of- 
ten resolved into knots (HH objects), moving with bulk velocities 
of 100-200 km s _1 and occasionally exhibit high degrees of colli- 
mation. The molecular outflows are slower, moving a few tens of 
km s _1 . The question of whether CO outflows originate from phys- 
ically distinct, wide-angle winds te.g. jKwan & TademanJl9 88) or 
whether they result from momentum transferred to molecular cloud 
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material from the jet (e.g., Konigl & Pudritz 2000) has not been an- 
swered iReip urth & Ba llvhoOltlLerv et all2002T) . 

Kwan & Tademari] <1988tll995l) review high spectral resolu- 
tion observations of the central regions of young stellar systems at 
the origin of collimated outflows. They point out that several char- 
acteristics of optical forbidden emission lines are most easily ex- 
plained by the presence of two winds of different origin. The two 
winds consist of the following: 1) a fast (> 100 km s _1 ) wind 
originating from the very central region of the accretion disc or 
from the star itself, and 2) a slow (~ 10 km s _1 ) wind launched 
from the outer (> 107?o) region of the accretion disc (the 'disc 
wind'). Evidence for at le ast two, physically distinct winds ex- 
ists for several objects (e.g. JSohj^FjohmT 993t lFridlund & Liseaul 
ll998UBacciotti et alfeOOOtlPvo et all2002l) . 

Prevalent theoretical models fo r winds launched from ac- 
cretion discs (see, e.g.. IShibata & Kudoblll999l iKonigl & Pudritzl 
2000, for a review) hold that the final wind velocity is of the or- 
der of the Keplerian rotational velocity of the launch point. The 
high velocity of observed optical jets suggests that they originate 
from a deep potential well iKwan & Tademarull98 8). requiring the 
launching region to be much less than an AU in extent, for reason- 
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able parameters. To within observational limits, the jets often ap- 
pear to emerge from their source with relatively large widths (sev- 
eral 10's of AU), but remain highly collimated immediately after- 
ward (e.g., Burrows et al. 1996; Eisloffel et al. 2000 ; Reipurfh et al. 
l2000t 120021) . We therefore adopt the view that jets are launched 
from a region within several stellar radii, initially with large open- 
ing angles (hereafter the 'central' or 'wide-angle' wind), but then 
become collimated within several 10's of AU. 

In wind launching theory, the self-collimation of the wide- 
angle wind (due to azimuthal magnetic fields c arried in the wind it- 
self) cannot occur for all of the wind material ( Blandf ord & Pavnel 
Il982t ISakurailll98^: IShu et aljll993) . In other words, without an 
external influence on the wind, the launching models predict that 
there is always some material leaving at very wide angles. On the 
other hand, theoretical studies of jets propagating into the surround- 
ing envi ronment most often assume that the flow is completely col- 
limated (Frank et al. 2002). Thus, there exists a gap in the theory 
between the launching and the propagation of YSO outflows. There 
is a need to understand how the initially wide-angle flows interact 
with their environment, presumably leading to a collimation of the 
entire flow. 

Recent work by iGardiner et alj 12002L also see references 
therein) considered the interaction of wide-angle winds with an 
unmagnetized, collapsing environment. They convincingly demon- 
strated that such flows may be collimated by shock focusing, espe- 
cially when azimuthal magnetic fields are present in the wide-angle 
wind. For the central wind, they adopted a flow similar to that pre- 
dicted by wind launching theory on large scales, though they were 
unable to directly model the launching of the flow. As a result, all of 
the ambient material, even that in the equatorial plane, was swept 
away in the wind, preventing any further accretion into the central 
region. This is in contrast to the expectation from a system con- 
taining an accretion disc, which prevents any outflow in a direction 
near the equatorial plane. In this paper, we consider the interaction 
of a wide-angle wind with a magnetized environment. Our method 
allows us to follow the details of the launching of the wide-angle 
flow from a region very near the central star, while still following 
the wind out to tens of AU, where it interacts with its environment. 

Since the central, wide-angle winds become collimated along 
the rotational axis of the accretion disc, the collimation process it- 
self must be associated with the angular momentum of the system. 
In this paper, we use numerical magnetohydrodynamic (MHD) 
simulations to explore the suggestion of lKwan & Tademar J ll988t) 
that disc-associated magnetic fields collimate a central, fast wind 
into an optical jet. Large-scal e disc fields may be generated 
by currents in the disc fepruit et alJ Il997l) . or they could be 
embedded in the surface of the disc and carried outward in a slow 
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initially vertical field to approximate the magnetic configuration 
expected from a disc on scales that are large compared to the wind 
launching region, but small compared to the size of the disc (e.g., 
Bland ford & Pavndll982L ISpruit et alJll997l) . The disc fields we 



consider are very weak compared to the ones that are likely to 
be involved in launching the central outflow (e.g., fields near the 
surface of the star). We will show that for typical YSO outflow 
parameters, disc fields of much less than 0. 1 Gauss are capable of 
completely collimating the flows. 

Our approach is t o use the 2.5D (a xisymmetric), resistive 
MHD code described by M att et ail 120021) to examine, parametri- 



cally, the effect of a vertical field on a wide-angle outflow. In order 
to gain insight into the physics involved, we first examine (in Sj2j 
the simplified situation of an isotropic, unmagnetized central wind 
interacting with a stagnant environment that contains a constant 
magnetic field. The main result is that collimation is inevitable, 
and that the size scale of collimation is inversely proportional to 
the strength of the field. In section|3] we assess the effect of a ver- 
tical field on a more realistic (anisotropic) outflow launched from 
the central region of a YSO. Our simulations are able to resolve the 
self-consistent interaction between the stellar dipole magnetic field 
and the inner edge of an accretion disc, which results in a highly 
structured, wide-angle flow from the central region. This central 
wind launching mechanism has been studied by several authors and 
is reviewed in section l3~T1 Using a system of nested computational 
grids, we follow this central flow to large distances, where it in- 
teracts with the vertical disc field. The interaction of this complex 
central wind with the disc field results in a highly structured flow, 
in contrast to the case with an isotropic central wind. We compare 
our results with observations, discuss the possible origin of a disc 
field, and suggest possible improvements to the model in section|4] 



2 COLLIMATION OF AN ISOTROPIC WIND 

In order to get a preliminary idea of the physics controlling the 
collimation of a central wind by a constant, vertical field, we first 
examine a simplified situation. Consider an isotropic wind that 
collides with a stagnant, constant density environment. In the ab- 
sence of any magnetic fields, the result will b e the w ell-studied 
case of a spherical wind -blown bubble (WBB; iKwok et allll978t 
iDvson & Williamsl H980). In that case, the wind acts as a piston 
pushing against the ambient medium, and an expanding structure 
results with two shocks and four distinct zones (free-flowing wind, 
shocked wind, shocked ambient, and undisturbed ambient mate- 
rial). 

The ambient material acts as an inertial barrier to the wind. 
When ordered magnetic fields are present, they introduce a direc- 
tional dependence to the inertial barrier. Consider, in particular, the 
simple case of a constant magnetic field permeating all of space 
and directed along the cylindrical z-axis (B z = constant). The wind 
flowing along the axis will experience no magnetic force in the z di- 
rection. However, the wind flowing perpendicular to the field (i.e., 
in the equatorial plane, z = 0) will be impeded by the magnetic 
pressure (assuming the wind is sufficiently ionized). If the field is 
relatively weak and unimportant near the source of the wind, it will 
be swept outward, in front of the wind. The dynamic pressure of the 
wind decreases with radius, so there will always exist a location at 
which the constant field becomes dynamically important and com- 
pletely stops the wind. In the case of a wind that is constant in 
time, continuity requires that the wind stagnating at the equator be 
deflected toward the z-axis. 

This relatively simple system, consisting of an isotropic cen- 
tral wind interacting with a plane-parallel mag n etic field, has 
been studied by a several authors iLevvl E22I lKonig]lll982t 
iFerriere et alll99ll:IStone & NormarJl992h . They were able to pre- 
dict the shape and expansion of the bubble of swept up ambient 
material, though they did not follow the details of the flow interior 
to that WBB. In this paper, we are interested in the collimation of 
the material originating in the wind, so this section contains a rela- 
tively simple study of this system and focuses on the conditions in 
the flow interior to the WBB. 

For simplicity, assume that the wind has a constant velocity 
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that is much higher than its own sound and magnetic wave (Alfven) 
speeds. Such is the case in most astrophysical outflows beyond sev- 
eral times the size of the launching region. Also, assume the exter- 
nal medium has a constant temperature and vertical magnetic field, 
B z , such that the Alfven speed is much higher than the local sound 
speed. Mathematically, 

Ww 2> C s ,w, Vw 3> VA,w, W,a 3> C s ,a, (1) 

where v v is the wind speed, c s is the sound speed, va is the Alfven 
speed, and the subscripts 'w' and 'a' refer to the value in the wind 
and ambient material, respectively. The conditions Q mean that 
the kinetic energy dominates in the free-flowing wind (thermal and 
magnetic energy are unimportant), and the magnetic energy domi- 
nates in the ambient medium (so that thermal pressure is unimpor- 
tant; also neglecting the effects of gravity, which is valid far from 
the wind acceleration region). 

Consider an initially stagnant medium and a spherical wind 
that turns on at t = 0. The wind will push against a weak magnetic 
field, 'carving out' a field-free region. The magnetic field will be- 
come important when the ram pressure in the wind (pv^, which 
falls of with distance from the star) is equal to the magnetic pres- 
sure of the external environment. For an isotropic wind, the balance 
of these two pressures on the equator (cylindrical z = 0) is equiva- 
lent to the balance of the ambient magnetic and wind kinetic energy 
densities. Therefore, the key parameter for determining the impor- 
tance of the magnetic barrier (B z ) to the wind is simply the ratio 

B 2 Z R 2 

M w «w 

where M v = 4TrR 2 pv v is the mass outflow rate in the wind. Note 
that a increases as the square of spherical radius (R) so that, even if 
B z is very small, magnetic forces will always be important at some 
distance. 

When the expanding shell of shocked material reaches a ra- 
dius at which a is near unity, the shell will deviate from spherical 
symmetry. The cylindrical r component of the wind velocity will be 
forced to zero. Since the wind travels faster than any information- 
carrying waves, the deceleration of the wind will be in the form of 
a shock. Pressure forces behind the shock will channel the post- 
shock flow toward the axis. The location of the deceleration shock 
in the equatorial plane is determined by setting a = 1 in Equation 
[2|and solving for R. Our simulations (presented in £|2,ll and H2.2t 
agree with this prediction for various wind and ambient parameters 
subject to conditions Q 



2.1 MHD Simulations 

To test the simple concept of an isotropic central wind interact- 
ing with a constant vertical field, and to develop insight into more 
complex problems, we have carried out a number of numerical sim- 
ulations. We use the 2.5D MHD code of lMatt et afl<2002l). and the 
reade r will find the details there (see also Goodson et al. 1997; Matt 
2002). The simulations employ a two-step Lax-Wendroff (finite- 
difference) scheme to solve the ideal MHD equations on a series of 
nested boxes. The equations allow for an ideal gas equation of state 
(7 = 5/3) and assume that the system is axisymmetric (d/dcj> = 
for all quantities). Each box consists of a 101 x 100 cylindrical 
(r-z) grid with constant grid spacing. The boxes are nested con- 
centrically, so that the inner box represents the smallest domain at 
the highest resolution. The next outer box represents a twice the 
domain size with half the spatial resolution (and so on for other 



boxes). For the simplified approach of this section, we use five 
boxes and ignore the effects of gravity and magnetic diffusion. Out- 
flow boundary conditions (all variables are continuous) exist on the 
right and top boundary of the outermost box, allowing material to 
leave the box without reflecting information inward. 

The wind source is a multiple layered boundary, consisting 
of a circle with an outer radius of 10.5 grid points in the smallest 
box, centered on the origin of the grid. The boundary conditions on 
the wind sphere are as follows: All quantities are fixed (velocity = 
0.0, B = 0.0, and density and pressure are constant) for R ^ 8.5 
(R = (r 2 + z 2 ) 1/2 ). For 8.5 < R < 9.5, the same is true except 
that B floats (i.e., it evolves naturally in the course of the simu- 
lations). For 9.5 < R ^ 10.5, the velocity is fixed at a constant, 
nonzero value and is directed radially, while all other quantities are 
the same as for 8.5 < R ^ 9.5. The existence of a layer with a zero 
velocity and a floating B interior to a layer with nonzero velocity 
is necessary to preserve V • B = near the win d source boundary 
(similar conditions are used bv lMatt et al . 2000). The wind source 
in the simulations represents a reference surface far from the re- 
gion where the wind is accelerated. In this way, the fact that the 
wind is supersonic and super-Alfvenic and the neglect of gravity is 
justified. 

The grid initially has a constant density, pressure, B z , and 
zero velocity. The key parameter is a*, where the asterisk sub- 
script refers to the value at the radius of our wind source, defined 
by cr* = B 2 J? 2 Mw '"w'- At t = 0, the wind begins to move out- 
ward from the wind source, pushing the ambient material ahead 
of it and sweeping up the ambient magnetic field. We first verified 
that, when B z = (a* = 0), the resulting WBB is spherically sym- 
metric. Nine other cases, with finite B z , were sufficient to cover a 
range of parameters. TableQcontains the parameter values of each 
case, and the simulation results are the topic of the next section. We 
consider the case with <r, = 0.0087, wind source to ambient density 
ratio (p*/p a ) = 10, wind Mach number (u w /c s , w ) = 5, and ambient 
Alfven to sound speed ratio (v A.a/c 8 . a ) = 5, as our baseline, with 
which to compare all other simulations. 



2.2 Results 

In the very earliest stage of the simulations, the interaction between 
the wind and the stagnant environment produces an expanding, 
spherically symmetric WBB. The flow has the following important 
features: 1) an interior filled with free-flowing wind, 2) a reverse 
shock (RS), where wind kinetic energy is converted into thermal 
energy, 3) a contact discontinuity (CD) separating the shock-heated 
wind gas from the shock-heated ambient gas, and 4) a forward 
shock where the expanding structure runs into ambient material and 
heats it. As the simulation evolves, and the structure expands, the 
vertical magnetic field (B z ) in the ambient material is swept aside 
in r. When the WBB reaches a location where a (Eqn.|2j is near 
unity, the RS slows and eventually stagnates. The position of the 
CD on the equatorial plane also stops, and pressure gradients in the 
region between the RS and CD channel material toward the axis. 
The end result is a steady-state, collimated flow (a jet) propagating 
along the z-axis. 

FigureQillustrates the flow for the baseline case near steady- 
state, within 40 wind source radii. The figure clearly demonstrates 
the existence of a RS at about 12 i?* (where the velocity of the 
wind suddenly drops) and the channeling via pressure gradients of 
the post-shock flow toward the axis. The deflected and compressed 
field lines (near the CD) provide a J x B force that balances the 



4 Matt, Winglee, & Bohm 



dynamic pressure of the wind. F igure [J] confirms the conceptual 
model of K wan & Tademanil ll988t compare with their fig. 2). 

We plot the equatorial stagnation location of the RS in seven 
of our simulations, as a function of the parameter a,, in figure [2] 
Square symbols represent the four cases with baseline ambient den- 
sity (top four rows in table0, though each include a unique value of 
B z . The triangles represent two cases with the ambient density de- 
creased by a factor of 10 (row 5 and 7 of table0, and the diamond 
is for the case with increased ambient density (row 6). We did not 
plot the results for the two simulations with various wind and ambi- 
ent temperatures (rows 8 and 9), since they were indistinguishable 
from the baseline case. The error bars in figure[2]represent the un- 
certainty due to the finite width of the shock, caused mainly by the 
relatively low numerical resolution. The dotted line represents the 
radial location where a = 1.0 for each a*. Figure [2] demonstrates 
that, within the uncertainty, the location of the RS on the equator is 
indeed predicted by equation|2| 

Figure [3] is a greyscale image of the pressure (logarithmic) 
for the baseline case. The data are from the same time as in figure 
Q but they represent the next outer box. The figure demonstrates 
the propagation of the jet along the axis. The jet is composed of 
shocked gas, and it is supported laterally by the balance of internal 
thermal pressure versus external magnetic forces. The outer edge 
of this jet (i.e., the edge of the high pressure region) represents the 
CD and is thus the location of the expanding bubble studied by 
(e.g.) IStone & Normanl ll992h . Though the jet propagates at a ve- 
locity that is higher than c Sja and comparable to «A,a, it is traveling 
slower than its own internal sound speed. In the more general case 
of an anisotropic central wind, it is possible for the shocked gas 
(the jet) to travel faster than its own internal sound speed. Also, if 
the gas were allowed to cool radiatively, the jet structure would be 
different. 

Since the flow is subsonic between the RS and CD, the pres- 
sure is roughly constant and equal to the decrease in the wind dy- 
namic pressure across the RS (i.e., the kinetic energy in front of 
the shock is converted to thermal energy behind it). This is evident 
in figure[3] The thermal pressure of the shocked wind extends the 
influence of the wind dynamic pressure to larger radii. The width 
of this extension on the equator (equal to the distance between the 
RS and the CD) depends upon how fast the shocked wind is moved 
toward the axis by small pressure gradients, which depends on how 
fast the material on the axis flows away (i.e., the speed of the jet). 

For given wind parameters and B z , the jet speed is determined 
by the conditions (in particular, the Alfven speed; see lLev\ll97lh in 
the ambient medium. The ambient density, p a , is important, since it 
behaves as an inertial barrier to the jet and it determines the Alfven 
speed. For a smaller p a , the jet moves faster, allowing the equatorial 
region between the RS and CD to be 'drained' more quickly by 
pressure gradients. Figure |4| shows the case that is identical to the 
baseline case but with an ambient density decreased by a factor of 
10 (row 5 of table0. The figure depicts an earlier time than figure 
Q because the jet moves more rapidly toward the top of the outer 
box (which limits the length of time the simulations can run), but 
a steady-state already exists within about 20 wind sphere radii. A 
comparison between figuresQand|4|reveals a dependence of the CD 
location (equivalent to where the magnetic field jumps in value) on 
the ambient density. 

The relationship between the CD location and p a is even more 
apparent in figure|5| a plot of the CD location versus the parameter 
a*. The symbols represent the same cases as in figure [2] and the 
dotted lines represent the radii of various values of a as a function 
of <t„ . For the cases of any given ambient density, the points follow 



the same trend (though shifted in position) as for the location of the 
RS. For an increased ambient density, the steady-state location of 
the CD on the equator moves outward. 



3 COLLIMATION OF A YSO OUTFLOW 

In this section we consider the more complex case of an anisotropic 
wind collimated by an initially vertical magnetic field embedded in 
an accretion disc. It is our goal to compare these results to observa- 
tions of outflows from YSO's, so we adopt a central wind launching 
mechanism that is appropriate for those systems. In the absence of 
any disc field, this MHD mechanism (reviewed in i|3.1t produces 
an outflow with a collimated (jet) component plus an uncollimated 
wind. Both components of the flow are relatively fast (100-200 km 
s _1 ) and of the order of the Keplerian rotation speed at the cen- 
tral region from which the wind is launched. Our simulations (de- 
scribed in £13.21 and results presented in i!3.3> resolve the detailed 
launching of the central wind, and they will show that a large-scale 
field, embedded in the accretion disc, will collimate the wide-angle 
component of the wind. This results in a more powerful and physi- 
cally broader jet than produced in the absence of the disc field. We 
also show that the simple results of the last section hold, even for 
the more complex outflow morphologies produced in this section. 

For simplicity, we will assume an initially constant, vertical 
field frozen-into the disc and permeating an environment above and 
below the disc that is initially in hydrostatic equilibrium. The small 
time-scales involved in properly capturing the launching of the cen- 
tral wind make it impractical for us to also follow the launching of 
disc plasma along the initially vertical field lines. Therefore, the ex- 
istence of a wide-angle, sl ° w i disc wind observed in some systems 
iKwan & TademanJl988T) is not addressed by our models. Such a 
wind may result from, or be responsible for, the existence of a large- 
scale, primarily vertical magnetic field (assuming it is sufficiently 
ionized; e.g.. lBlandford & Pavndll982h . Our simple treatment of 
the disc field will, in some ways, emulate a magnetized disc wind 
that can act as a chan nel for the interior, wide-angle wind (e.g., 
lOuved & Pudritzll 19971) . We propose that the relatively broad jet, 
produced by the collimation of a central wind by this disc field, can 
explain some of the observational properties (e.g., the mass flux, 
velocity, energy and momentum budget, and collimation) of opti- 
cal jets. 



3.1 The Central Wind 

For th e central wind, we ado pt a model stu died by H ayashi et alJ 
h99dl . lGoodson et aljjl997l) . and Miller & StoneHl997l) . in which 
a rotating star couples to the inner edge of a conductive accretion 
disc via a stellar dipole magnetic field. In this model, differential ro- 
tation between the star and disc twists up the magnetic field, adding 
a component. The magnetic pressure associated with causes 
the stellar magnetosphere to expand (or 'inflate') above and below 
the disc. Plasma loaded onto the inflating field lines is driven out- 
ward in a burst. When the inflation speed becomes larger than any 
plasma wave speeds, the field lines become effectively open, and 
material is magnetocentrifugally launc hed from the disc inn er edge 
via the Blandford-Payne mechanism I Blandford & Payne 1982). 
Plasma launched from the disc inner edge carries away angular mo- 
mentum, so the disc spins down and moves inward. The inward mo- 
tion of the disc forces together the open field lines from the star and 
disc of opposite polarity, instigating a reconnection of those lines. 
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After reconnection, material at the disc inner edge is again mag- 
netically coupled to the star, but with less angular momentum than 
before, so that it accretes via funnel flow along field lines onto the 
star. The accretion of material unloads the stellar magnetosphere, 
allowing it to expand outward and diffuse into the new inner edge 
of the disc. When the field begins to thread the disc inner edge, dif- 
ferential rotation again twists up the field, and the process repeats. 
This process of episodic magnetospheric inflation (EMI) regulates 
both the accretion and ejection of material within the region of the 
disc inner edge. 

B y following the outflow to several AU. lGoodson et alii 19971 
1999) showed that the EMI mechanism produces a jet (hereafter, 
the 'EMI jet') that has knots associated with each reconnection 
event. Note that the typical oscillation period for the EMI mech- 
anism is about 20-30 days, t oo short to explain the knot spacing in 
observed jets (e.g., in HH 30 lBurrows et alll996h . The rest of the 
flow travels out at a wide opening angle, shocking the surface of 
a flared accretion disc. The flow speed for the EMI jet and wide- 
angle wind compone nt is > 100 km s~ , in rough agreement with 
observations ( Reipurt h & Ballvl200ll) . The energetics and the mass 
outflow rate in the EMI wind are also within the range observed in 
jets, but most of the kinetic energy and mass in the EMI wind is 
contained in the wide-angle component (not in the jet). In addition, 
the width of the EMI jet compone nt is a few AU, an or der of mag- 
nitude smaller than observed (e.g., B urrows et all 19961) . 

In their lim ited parameter study of the EMI mechanism, 
iMatt et all (2002) showed that a weak magnetic field, aligned with 
the rotation axis and threading the accretion disc, did not affect 
the launching of the EMI outflow. This should be true, as long as 
the stellar dipole field is stronger than the vertical field at the in- 
ner edge of the disc. Their work focused on the integrated outflow 
properties, rather than the details of the wind morphology, and they 
were unable to address the collimation of the outflow by the vertical 
fields because the size of their largest simulation grid reached only 
to 0.7 5 AU. In this section, we follow up the work of Ma tt et alJ 
(2002) and, using insight gained from section|2| we will show that 
when one follows the outflow to large enough distances, even a 
weak field will collimate the enti re EMI outflow. 

Many other authors (e.g.. Ish u et alJ 1 19941 lLovelace et alJ 
ll995tlFendt & ElstneJl999tlKuker et all2003l) have discussed out- 
flows produced by the interaction between a stellar magnetosphere 
and the inner edge of the disc. The EMI mechanism operates only 
when the magnetic diffusivity is relatively small in these systems 
iGoodson & WingledlT999l) . but the outflows produced in all of 
these studies contains a component that flows out at wide angles. 
We choose to work with the EMI model as one example, but our 
conclusions apply to any wide-angle, central flow. 

3.2 Simulation Method 

As in section^ we use the 2.5D MHD code of IMatt et all <2002h . 
For the more complex simulations carried out in this section, the 
code solves the ideal MHD equations with the added physics (via 
source terms in the momentum, energy, and induction equations) 
of gravity and Ohmic diffusion. We also adopt their baseline sim- 
ulation case, in which a 1 Mq, 1.5 Rq T-Tauri star with a 1.8 day 
rotation period interacts with a flared accretion disc (after the a- 
disc model of Shakura & Sunvaev 1973) via an axis aligned dipole 
magnetic field (= 909 Gauss on equator). The reader will find more 
details of the simulations and initial conditions in IMatt et alji2002T) 
and Goodson et al. ( 1997, 1999). Our simulations differ from those 
of Ma tt et alJ 120021) only in that we use ten nested boxes, while 



they used only four. Each simulation box consists of a 401 x 100 
cylindrical (r-z) grid, and the boxes are nested concentrically (see 
EH . In this way, we fully capture the launching of the outflow via 
the EMI mechanism in the innermost box, while our outermost box 
follows the flow out to 24 AU. 

First, we simulated the baseline case and follow the outflow to 
large radii, allowable by the increased number of boxes. We then 
ran four other cases with an additional vertical and constant mag- 
netic field (B z ) initialized everywhere in the simulation region. All 
of the cases are identical except for the strength of the initial ver- 
tical field. So the baseline case has B z = 0, and the others have 
\B Z \ = 0.050, or 0.025, Gauss. The vertical field can interact (e.g., 
reconnect) with the field carried in the outflow, so there is a need 
to assess the effect of polarity of the vertical field. Thus, for each 
case with nonzero \B Z \, we ran two cases, each with opposite po- 
larity. Table |2| lists the values of B z , the simulation time at which 
the simulations were stopped (if), and some notes for each case. 

Our choice of vertical field strengths was based on the need for 
the field to be strong enough to collimate the flow within the outer 
box and weak enough not to effect the EMI mechanism (which 
gives an upper limit of ~ 0.1 G). We used equation|2|to estimate 
the required field. 

3.3 Results 

In the absence of any disc field, the basic properties of the EMI 
wind include 1) a relatively high density, narrow jet (the EMI jet) 
along the z-axis and 2) a wide-angle wind. These features are ap- 
parent in Figure [6] a snapshot from the tenth simulation box of 
the outflow in the B z = baseline case at 241 days. The star is at 
r = z = 0, and the data has been reflected about the r = axis 
to better demonstrate the flow. The Figure also contains data from 
the eighth and ninth simulation boxes, which extend to 6 and 12 
AU, respectively, in z. The outflow is highly structured, contain- 
ing knots in the jet and 'wispy' structures in the wide-angle wind. 
These structures are correlated with oscillations of the EMI mech- 
anism. The entire flow originates from a region of < 30Rq , where 
the inner edge of the disc is located. By the time of this snapshot, 
material originating in the disc or on the star has traveled to ~ 16- 
18 AU, and the flow at larger z is transient (i.e., due to the initial 
expansion into a stagnant environment). 

FigureQis a snapshot from the B z — 0.025 G (left panel) and 
B z — —0.025 G (right panel) cases at 241 and 232 days, respec- 
tively. It is evident that the wide-angle flow has been completely 
collimated by the disc-associated field to a width of 6-7 AU. In 
other words, the EMI jet and wide-angle wind have been combined 
into a broad, collimated structure (hearafter, the 'broad jet'). Note 
that the narrow EMI jet is still apparent along the axis, and a rela- 
tively dense shell has formed on the outer edge (in r) of the broad 
jet, where material traveling at large polar angles piles up against 
the disc field before changing direction. However, the central EMI 
jet and the outer shell exists only to a height that is comparable to 
the collimation height of the broad jet (at z ~ 12 AU). A compari- 
son between the left and right panels of Figure Qreveals that there 
is almost no difference in collimation radius for the cases with op- 
posite disc field polarity. As in figure|6| the main flow has reached 
to z ~ 16-18 AU. Since this distance is comparable to the width 
of broad jet, the final collimation width may become a little larger 
than achieved by the time of this snapshot. 

Figure|8|is a snapshot from the B z = 0.05 G (left panel) and 
B z = —0.05 G (right panel) cases at 244 and 196 days, respec- 
tively. This figure shows the collimation of the entire central flow 
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to a width of 3-5 AU. Again, the two panels show that there is al- 
most no difference in collimation radius between the two cases of 
opposite disc field polarity. As in figure^ both panels also contain 
the central EMI jet and a dense shell enclosing the broad jet to a 
height of z w 8 AU, where the broad jet becomes collimated. 

The data plotted in figuresQand[8]are consistent with the pre- 
diction of equation |2| that the collimation radius is inversely pro- 
portional to the disc field strength. For a more quantitative analysis, 
we plot various outflow properties as a function of polar angle (9) 
at a constant spherical radius in figuresl9 hl II Because the wind is 
episodic and highly structured, and because our time sampling for 
data outputs is fairly course (2.8 days), it is necessary to represent 
the outflow properties as a time-averaged quantities that covers sev- 
eral periods of the EMI oscillation (which is around 20-30 days). 
We therefore calculate the outflow quantities for each case at 90 
gridpoints out in our eighth simulation box (= 5.4 AU) and average 
them from 112 days until the end of the simulation. 

Figure|9^ shows the mass flux, $m{8), for all cases. The ex- 
istence of the narrow EMI jet at 8 < 10 degrees is represented 
by the strongly peaked mass flux there for all cases. Note that, for 
8 > 65°, $a/ becomes negative, as this represents material inside 
the accretion disc. Figure |9jj shows the mass outflow rate per dif- 
ferential polar angle; dM/dd = 2ttR 2 &m(8) sin#, where R is 
the spherical radius. It is evident that, while $m is greatest near 
the pole, most of the mass in the flow exists at large opening an- 
gles (8 > 25°). For every case, there is a peak in <I>ai and dM/dd 
where 8 is between 25° and 55°. For larger B z , this peak is de- 
flected to smaller angles. The wide-angle flow is not completely 
collimated at z = 5.4 AU for any of the cases (see figs. [7] and [8}, 
but the deflection of this peak indicates the onset of collimation. 
The cases with positive B z are shifted a few degrees less than their 
negative B z counterparts. A positive and negative B z corresponds 
to the stellar magnetosphere havi ng an initially op en and initially 
closed topology, respectively (see lMatt et aljEo02l) . This suggests 
that the open topology is slightly less efficient at collimating at 5.4 
AU (presumably due to the release of magnetic energy from mag- 
netic reconnection between disc field lines and those carried in the 
outflow), though the final collimation radius is not much affected 
by the topology of the field (as evident in figs. Q and Finally, 
there is a bump in $m between 55° and 65° for the cases with 
\B Z \ > 0. This feature is much less prominent in the momentum 
and kinetic energy (see figs. HOI and llH of the outflow, indicating 
that it has relatively high density and low velocity compared to the 
rest of the outflow. It represents material being pushed aside by the 
compression of disc field lines. 

The total mass outflow rate can be obtained by integrating 
dM/d8 over a desired range of 8. For all cases, the positive mass 
outflow rate (excluding the bump at 55° < 8 < 65° for the 
\B Z \ > cases) is ~ 1.5-2.0X 1O _8 M0 yr _1 in one hemisphere. 
This agrees very well with the typical val ue of M ~ 3 x 10~ 8 Mp 
yr _1 derived from observations of jets iReipurth & Ballvjl200ll> . 
Note, however, that for the B z = case, the mass loss happens to- 
ward large opening angles, while for the \B Z \ > cases, the whole 
of M is contained in the broad jet. 

All cases also show a very strongly peaked EMI jet in the 
momentum flux (Fig. llOh ). p(v 2 + v 2 ) x ^ 2 (t; • R), and kinetic en- 
ergy flux (Fig.fTob). 0.5p(v 2 + v'i)(v ■ R). The total kinetic en- 
ergy outflow rate (mechanical luminosity) of the flow is ~ 10 32 
ergs s~ , lower t han th e typical value of ~ ILq reported by 
iReipurth & BallvlfeOOll) . The wide-angle peak in both fluxes of fig- 
ure ^)| follow the same case-by-case behavior as in figure [9] That 
is, the peak shifts to smaller angles for stronger \B Z \. 



Figure \H] shows the linear momentum outflow rate (thrust) 
per differential polar angle; dp x /dd = 2tvR 2 s'md x p x -flux(8), 
where p x is m (v 2 +vj) 1/2 forfttt. mv r forfTTt. and mv z forfTTl. 
The total vertical momentum outflow rate lilt ') is ~ 10 25 dynes 
for all case s, lower than the typical value of ~ 10 26 gm cm s~ 2 , 
reported by Reiourth & Ballvl 1200 ll) . In all cases, the total thrust 
il H i) is dominated by the vertical thrust till :) for 8 < 40°. For 
the baseline case, the r-directed <1 lb ) and vertical thrust 1111 :) are 
comparable for 8 > 40°, indicating that the flow is moving in a 
direction of ~ 45° there. This is because the wide-angle component 
of the B z — case is not collimated. However. I lib demonstrates 
that, for increasing \B Z \, the total r-directed thrust decreases. This 
is because the expanding wind experiences a negative-r-directed 
J x B force from the perturbed disc field. The peak in vertical 
thrust (c) shifts to smaller polar angles for larger j B z j (indicating 
the onset of collimation), but the total vertical thrust (integrated 
over 8) is roughly the same for all cases. 

The calculation of the outflow thrust (Fig. 1111 allows us to 
apply the results of section|2|quantitatively, though roughly, to the 
EMI outflow. Since it is only the r-directed component of the wind 
that deflects the disc field before the collimation radius is reached, 
we will substitute M w Vw of equation|2|with p r and set a equal to 
unity. The collimation radius can then be estimated as 



The p r (integrated over 8) for the B z = case is ~ 6.5 x 10 
dynes in one hemisphere. This predicts a r co i of 4.8 and 9.6 AU for 
the \B Z \ — 0.05 and 0.025 G cases, respectively. This prediction 
is slightly larger than the simulated value for the \B Z \ — 0.025 
cases (see figs. Q and table |2j, but remember that the simulated 
flow in that case has not quite reached its full collimation radius. 
Also, for a magnetized central wind (as in the EMI outflow), there 
will be an increase in the collimation (decrease in r co i) as B^, car- 
ried i n the wind is co mpressed in the outer shell (see figs.0and|8j 
iGardiner et alj2002t) . 



4 SUMMARY AND DISCUSSION 

The relatively high velocity of observed YSO jets requires that jet 
material is launched from deep within the gravitational potential 
well (i.e., within several stellar radii). In order to explain the ob- 
served jet widths and their degree of collimation, the jet material 
must be launched initially with wide opening angle and become al- 
most completely collimated at a distance from the star that is com- 
parable to the width of the jet. Self-collimation of the wide-angle 
wind (due to magnetic fields carried in the wind itself) cannot oc- 
cur for 100% of the wind material; there is always some material 
leaving at very wide angles. This argues for the existence of a col- 
limator that is external to the central wind. We have explored one 
such possibility. 

Specifically, we used time-dependent MHD simulations to ex- 
tend the work of iKwan & Tademarul {l988) and demonstrate that 
any outflow (even if initially isotropic) will be eventually colli- 
mated by a constant, parallel magnetic field of any strength. The 
collimation of the central wind is inevitable, assuming that the wind 
is sufficiently ionized to interact with the field. This is because the 
ram pressure in the wind always decreases with radius for a diverg- 
ing wind, and the magnetic pressure is constant. The collimation 
radius is roughly given by the location where the two energy den- 
sities balance. 
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In section [5] we applied this concept to realistic YSO out- 
flows by first adopting the episo d ic mag netospheric inflation 
(EMI) model faavashi et alJ Il994 lOoodson et alJ Il997t 11994 
iMiller & Stone 1997) for launching the outflows. We then assumed 
that the entire simulation region was initially threaded by a vertical 
magnetic field. The result is that the wide-angle component of the 
EMI wind is collimated by disc fields of <C 0.1 Gauss. In this way, 
the entire mass, momentum, and energy budget of the central out- 
flow is channeled into a relatively fast (100-200 km s _1 ), broad jet 
with a width of several AU. The channeling occurs within a vertical 
height of the order of the jet width, consistent with observational 
constraints (see iQ. 

The morphology of the resulting broad jet includes a relatively 
narrow density enhancement along the axis. This feature is the nar- 
row jet produced by the EMI mechanism alone (i.e., even when 
there is no disc field). The collimation of the wide-angle wind re- 
sults in a relatively high density shell of shocked, outflowing ma- 
terial surrounding the flow out to at least the height of collimation. 
The EMI jet and shocked shell together resemble a 'pitchfork' in 
proj ection (see fi gs . and |8j ■ 

iBallv et all 120 03) recently suggested that x-ray emission from 
a region C/.'5 to l'.'O west-southwest of the YSO L1551 IRS 5 may 
be explained by shocks associated with the collimation of the main 
outflow from that region. Using kinematic information from IR slit 
observations of that outflow, Pvo et al. 1 2002) demonstrated the ex- 
istence of two different flows (not including the molecular outflow 
discovered bv lSnell et alJll9 80): an ionized, highly collimated jet 
moving at ~ 440 km s _1 , and a partially ionized, wide-angle wind 
moving at ~ 200 km s . The wide-angle wind shows a kine- 
matic sign ature of inc reased collimation as it moves further from 
the source. |pvo et alJ (2002) also noted a pitchfork-like morphol- 
ogy in an I band image of the outflow, which they suggested could 
represent the limb-brightened edges of the wide-angle wind under 
collimation. Note, however, that the southernmost of the pitchfork 
'tines' has been interpreted as a separate jet from a companion star 
iFridlund & Liseaull998tlrlartigan et alj2000h . 

In this paper, we have assumed the special geometry of an 
axis-aligned, vertical field, which requires that the field is as- 
sociated with the accretion disc. This field could represent a 
field carried in a disc wind (moving more slowly than the cen- 
tral wind), though we don't follow such a wind in our simu- 
lations. Specifically, if the disc is magnetized, and if there is 
a magnetized outflow from the disc, the field will roughly fol- 
low the wind streamlines in the poloidal plane. Indeed, the pre- 
dicted p oloidal fi eld geo metry f or disc wind launching mod- 
els (e.g.. lUchida & Shibat3ll985l: IShibata & Uchidj[i986l Il99fi 
lOuved & Pu^t^^9^nUstvu^ovaet^L^^99l) is essentially verti- 
cal near the rotation axis. 

Disc wind models also predict that the field has a helical struc- 
ture, due to the winding up of field lines by the rotation of the disc. 
The azimuthal field, added in this way, increases the co llimation of 
the disc wind telandford & Pavnell982HSakuraill98^) . and would 
therefore further constrain any central, wide-angle flow that inter- 
acts with it. Due to the long time-scale for the winding of disc field 
compared to the time-scale for the launching and propagation of 
the central wind by the EMI mechanism, our simulations do not 
capture the effect of the azimuthal component expected in the disc- 
associated field. 

Also, due to differential rotation and possibly turbulence in 
a real accretion disc, it is reasonable to assume that any field 
present in the disc may b e disordere d, and possibl y chaotic (e.g., 
iBlandford & Pavndll982t lHawlev & Balbuslll992t) . A disc wind 



would therefore be threaded by a magnetic field with direction re- 
versals at irregular intervals (i.e., flux tubes with opposite polarity). 
Such direction reversals do n ot affect the qualitative behavior of the 
field (see, e.g., Tsinganos & Boaovalov 2000). That is, whether the 
vertical field has a constant or chaotic polarity, it will always act 
to collimate the flow (see also Li 2002 for similar description with 
chaotic field). 



The interaction of a wide-angle wind with an external, disc- 
associated magnetic field (as presented here) is not the only possi- 
ble larg e-scale co llimation mechanism for YSO outfl ows. A few 



authors (Frank_&Noriega-Crespo 


19941: iFrank & Mellema 


1996 


iLi & Shdll996t iMellema & Frank 


19971 iDelamarter et al 


2000 


Gardiner et all20021) have discussed the interaction of wide-angle 



flows with an unmagnetized environment (see also iLerv et alJ 
2002). They have found that the wide winds may be completely 
collimated by shock focusing, especially when azimuthal magnetic 
fields are present in the wide-angle wind. It is not clear what the ob- 
servational signature of such an interaction would be (for possible 
detections, see lMatt & B6hmll2003l:lBallv et all2003l) . but it takes 
place relatively close to the source and therefore may be difficult to 
disentangle from the complex emitting environment there. 
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Table 1. Isotropic wind simulation parameters. 
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Table 2. YSO simulation parameters. 



B z 


tf 


Notes: 


(Gauss) 


(days) 




0.000 


241 


no disc field 
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241 


r'col > 6-7 AU 


-0.025 


232 


J-col > 6-7 AU 


0.050 


244 


r co i ?s 3-5 AU 


-0.050 


196 


r C ol ~ 3-5 AU 




r (R«) 

Figure 1. Greyscale image of log density (darker is higher density), field 
lines, and velocity vectors illustrate the near steady-state flow for the <r» = 
0.0087, baseline case. 
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Figure 6. Grayscale image of density (logarithmic) and velocity vectors from the case with no vertical field at 241 days. Log p > —15.3 gm cm 3 is black, 
log p < —19.3 gm cm~ 3 is white. 
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Figure 7. Greyscale images of density (logarithmic), magnetic field lines, and velocity vectors from the case with a 0.025 Gauss vertical field. The two panels 
represent simulations that are identical except for the polarity of the vertical field (left panel has positive B z , the right has negative, and the time in days is 
shown in each panel). The greyscale is identical to figure^] and the line spacing is proportional to field strength. 
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Figure 8. Greyscale images of density (logarithmic), magnetic field lines, and velocity vectors from the case with a 0.05 Gauss vertical field. The two panels 
represent simulations that are identical except for the polarity of the vertical field (left panel has positive B z , the right has negative, and the time in days is 
shown in each panel). The greyscale is identical to figure|5] and the line spacing is proportional to field strength. 




Figure 3. Greyscale image of log pressure (darker is higher pressure), field 
lines, and velocity vectors illustrate the flow at a late time and a larger scale 
for the a, = 0.0087, baseline case shown in figurelTI 
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Figure 4. Greyscale image of log density (same greyscale as fig.QJ, field 
lines, and velocity vectors illustrate the flow at near steady-state for the 
<r* = 0.0087 case with ambient density a factor of 10 less than the baseline 
case (fig.fT1andl3l. 
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Figure 9. a) Time-averaged mass flux and b) time-averaged mass outflow 
rate per radian (dM / d8, see text) versus polar angle, at a constant radius of 
5.4 AU. The solid line in each panel represents the case without a vertical 
field (B z = 0), the dotted lines are for cases with \B Z | = 0.025 G, and the 
dashed lines are for cases with \B Z | = 0.05 G. Bold lines are for B z < 0, 
and the lighter lines are for B z > 0. 



Figure 5. Equatorial location of the jump in B z as a function of cr, . The 
simulated data (symbols for various cases indicated on figure) follow the 
predicted trend for the standing shock (dotted lines for different values of 
a, chosen to fit the data) but also reveal a dependence on conditions of the 
ambient material. The errors in the simulated values (vertical bars) are due 
to the finite width of the jump in B z (due mainly to numerical resolution). 
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Figure 11. Time-averaged momentum outflow rate (thrust) per radian 
(dp /d9, see text) for a) the total, poloidal momentum, b) the radial compo- 
nent, and c) the vertical component versus polar angle, at a constant radius 
of 5.4 AU. All three panels are plotted on the same scale. The line styles 
represent the same cases as in figurel9l 



